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Abstract 

The shapes of jets with transverse energies, Ej^^, up to 45 GeV produced in neutral- 
and charged-current deep inelastic e+p scattering (DIS) at > 100 GeV^ have been 
measured with the ZEUS detector at HERA. Jets are identified using a cone algorithm 
in the i] — Lp plane with a cone radius of one unit. The jets become narrower as E^r^^ 
increases. The jet shapes in neutral- and charged-current DIS are found to be very 
similar. The jets in neutral-current DIS are narrower than those in resolved processes 
in photoproduction and closer to those in direct-photon processes for the same ranges in 



Ej:^ and jet pseudorapidity. The jet shapes in DIS are observed to be similar to those in 



e^e~ interactions and narrower than those in collisions for comparable Ej^^. Since the 
jets in e^e~ interactions and e"*"p DIS are predominantly quark initiated in both cases, 
the similarity in the jet shapes indicates that the pattern of QCD radiation within a quark 
jet is to a large extent independent of the hard scattering process in these reactions. 
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1 Introduction 



The internal structure of a jet is expected to depend mainly on the type of primary parton, 
quark or gluon, from which it originated and to a lesser extent on the particular hard scatter- 
ing process. For cone jet algorithms 0, 0] a useful representation of the jet's internal structure 
is given by the jet shape ||^. At sufficiently high jet energy, where fragmentation effects be- 
come negligible, the jet shape should be calculable by perturbative Quantum Chromo dynamics 
(pQCD). pQCD predicts gluon jets to be broader than quark jets as a consequence of the gluon- 
gluon coupling strength being larger than that of the quark-gluon coupling [Q. Measurements 
of the jet width in e^e~ interactions at LEPl have shown that gluon jets are indeed broader 
than quark jets 0. The dependence of the structure of quark and gluon jets on the production 
process can be investigated by comparing measurements of the jet shape in different reactions 
in which the final-state jets are predominantly quark or gluon initiated. 

Measurements of the integrated jet shape were made in pp collisions at a/s = 1.8 TeV 
using charged particles as well as both neutral and charged particles 0], and a qualitative 
agreement with 0(a|) QCD calculations P, |^ was found. Measurements of the integrated and 
differential jet shapes were made in e^e~ interactions at LEPl using both neutral and charged 
particles and were found to be well described by leading-logarithm parton-shower Monte 
Carlo calculations. It was observed |^ that the jets in e^e~ are significantly narrower than 
those in pp and most of this difference was ascribed to the different mixtures of quark and gluon 
jets in the two production processes. 

Measurements of the integrated jet shape in quasi-real photon proton collisions at HERA 
have recently been presented [0 and were found to be well described by leading-logarithm 
parton-shower Monte Carlo calculations except for the inclusive production of jets with high 
jet pseudorapidity (77-'^*) and low jet transverse energy {Ej^^). Fixed-order perturbative QCD 



calculations at the parton level |TT[ are able to describe the measured jet shapes within the 
uncertainties on the matching between the theoretical and experimental jet algorithms. 

At HERA, jet production has been observed in both neutral- [l^, |13| and charged-current 



H| deep inelastic ep scattering (DIS) at large (where is the virtuality of the exchanged 



boson). In this paper, measurements of the differential and integrated jet shapes in neutral- and 
charged-current DIS at > 100 GeV^ are presented. The data sample used in this analysis has 
been collected with the ZEUS detector in e~^p interactions at the HERA collider. To compare 
with measurements of the jet shapes in pp, 'jp and e~^e~ collisions, jets are searched for with an 
iterative cone algorithm [10| with radius i? = 1 in the pseudorapidityQ (77) - azimuth (y?) plane 
of the laboratory frame. Jets have been selected with jet transverse (with respect to the proton 
beam direction) energy E'^f* > 14 GeV and jet pseudorapidity in the range — 1 < rj^'^^ < 2. The 
jet shape has been measured using the ZEUS calorimeter and corrected to the hadron level. 
The measurements are presented as functions of Eij^^ and 77-^^*. The measured jet shapes are 
compared to similar measurements in other reactions and to leading-logarithm parton-shower 
Monte Carlo calculations. 



^The ZEUS coordinate system is defined as right-handed with the Z-axis pointing in the proton beam 
direction, hereafter referred to as forward, and the X-axis horizontal, pointing towards the centre of HERA. 
The pseudorapidity is defined as 77 = — ln(tan | ) , where the polar angle 9 is taken with respect to the proton 
beam direction. 
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2 Experimental setup 



During 1995 and 1996 HERA operated with protons of energy Ep = 820 GeV and positrons of 
energy = 27.5 GeV. The ZEUS detector is described in detail in The main subdetec- 

tors used in the present analysis are the central tracking system positioned in a 1.43 T solenoidal 
magnetic field and the uranium-scintillator sampling calorimeter (CAL). The tracking system 
was used to establish an interaction vertex and to select neutral- and charged- current DIS 
events. The CAL is hermetic and consists of 5918 cells each read out by two photomultipliers 
tubes. Under test beam conditions the CAL has energy resolutions of 18%/\/E for electrons 
and 35% /^/E for hadrons. Energy deposits in the CAL were used to identify the scattered 
positron, to find jets and to measure jet energies. Jet energies are corrected for the energy 
lost in inactive material in front of the CAL. This material is typically about one radiation 
length. The effects of uranium noise were minimised by discarding cells in the electromagnetic 
(EMC) or hadronic (HAC) sections if they had energy deposits of less than 60 MeV or 110 MeV, 
respectively. A three-level trigger was used to select events online W^, IB, O . 



3 Data selection 

Neutral-current (NC) DIS events have been selected offline from the ZEUS 1995 data sample, 
which corresponds to an integrated luminosity of 6.3 pb~^, using criteria similar to those re- 
ported previously |T^, The main steps are briefly discussed here. The scattered positron 



candidate has been identifled by using the pattern of energy deposits in the CAL [T^. The 



energy (-Eg') and polar angle {9^') of the positron candidate have been determined from the 
CAL measurements. The variable has been reconstructed by the double-angle method 
(Q'da) [13] 5 which uses 9e' and an angle that corresponds to the direction of the scattered quark 
in quark-parton model type events. This second angle has been determined from the CAL 
measurements of the hadronic flnal state. The following requirements have been imposed: 

• A positron candidate of uncorrected energy E(.> > 10 GeV. This cut ensures a high and well 
understood positron flnding efflciency and suppresses background from photoproduction 
events, where the scattered positron escapes down the rear beampipe. 

• He < 0.95, where j/e = 1 — Ef.i{l — cos6e')/{2Ef,). This condition removes events where 
fake positron candidates are found in the forward region of the CAL. 

• The total energy not associated with the positron candidate within a cone of radius 0.7 
units in the f] — <f plane around the positron direction must be less than 5 GeV. This 
condition removes photoproduction and DIS events where part of a jet has been falsely 
identifled as the scattered positron. 

• A track is required to match the positron candidate identifled in the CAL for r]e < 2, 
where rje is the pseudorapidity of the positron candidate. This requirement suppresses 
cosmic rays, beam-halo muons, photoproduction and DIS events where an electromagnetic 
shower in the CAL has been falsely identifled as the scattered positron. 

• For rjf, > 2 the transverse energy of the positron candidate should be larger than 20 GeV. 
This requirement further reduces the number of fake positrons in the forward region of 
the CAL. 
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38 GeV < {E — pz) < 65 GeV, where E is the total energy as measured by the CAL, 
E = Y,i Ei, and pz is the Z-component of the vector p = J2i EiVi] in both cases the sum 
runs over all CAL cells, Ei is the energy of the calorimeter cell i and rj is a unit vector 
along the line joining the reconstructed vertex and the geometric centre of the cell i. This 
cut removes events with large initial-state radiation and further reduces the background 
from photoproduction. 

Events have been removed from the sample if there was a second positron candidate 
with energy above 10 GeV and without a track match, and the energy in the CAL 
after subtracting that of the two positron candidates is below 10 GeV. This requirement 
removes elastic Compton scattering events (ep Gjp). 

ft/VE't < 3 GeV^/^ „here ft is the missing transverse momentum as measured with 



the CAL (ft = \Ipx +Py) ^"^^ Et is the total transverse energy in the CAL. This cut 
removes cosmic rays and beam-related background. 

• The vertex position along the beam axis must be in the range —30 < Z < 36 cm. 

• Ql^ > 100 GeV2. 

Charged-current (CC) DIS events have been selected offline from the ZEUS 1995 and 1996 
data samples, which correspond to an integrated luminosity of 14.8 pb~^, using criteria similar 
to those reported in [Q. The variable has been determined using the method of Jacquet- 
Blondel (Q'jb) IMHj which uses the information from the hadronic energy flow of the event. The 



following conditions have been imposed: 

• > 11 GeV. This cut ensures high trigger efficiency. 

• ft/ Et > 0.5. This cut rejects photoproduction and beam-related background. 

• The vertex position along the beam axis should lie in the range —30 < Z < 36 cm. 

• At least one track should point to the vertex. This requirement rejects cosmic rays and 
beam-gas interactions. 

• The number of tracks not associated to the vertex must be less than 20% of the total 
number of tracks. This cut further reduces the background from beam-gas interactions. 

• The difference A(p between the azimuths of the net transverse momentum as measured by 
the tracks associated with the vertex and as measured by the CAL has been required to 
fulfill I Ay? I < 1 rad. This requirement removes overlays of cosmic rays on ep interactions. 

• piracks^^^ > 0.1, where P*"^"-^^^ is the net transverse momentum of the tracks associated 
with the vertex (this condition has not been applied if f > 25 GeV). This cut rejects 
beam-related background, in which f is pointing to small polar angles, and events with 
additional non-ep related energy deposits in the CAL (mainly cosmic rays). 

• The event has been removed from the sample if there was an isolated positron candidate 
with energy above 10 GeV and a track match. This condition removes NC DIS events. 

• Pattern recognition algorithms based on the topology of the CAL energy distribution 
were applied to reject cosmic rays and beam-halo muons. 
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• Qlon- > GeV^, where Q^orr denotes the corrected value of Qj^ as described in |T^. 
The resolution in the reconstruction of is ~ 25%. 

A search for jet structure using the CAL cells (see Section 5) has been performed on both 
samples (NC and CC DIS), and events with at least one jet of 'corrected' transverse energy 
(see Section 5) E^* > 14 GeV and —1 < 77-''^* < 2 have been retained. The selected sample 
of NC (CC) DIS consists of 6926 (231) events containing 7092 (233) jets. In both cases, the 
background from photoproduction has been estimated using Monte Carlo techniques and was 
found to be below 1%. 

4 Monte Carlo simulation 

The response of the detector to jets and the correction factors for the jet shapes have been 
determined from samples of Monte Carlo (MC) events. 

NC and CC DIS events have been generated using the lepto program |2^ interfaced to 



HERACLES via DJANGO Hj. The HERACLES program includes photon and Z" exchanges 



and first-order electroweak radiative corrections. The CTEQ4D |^5[ NLO proton parton den 
sities have been used. The hadronic final state is simulated using the colour-dipole model 



including the leading-order (LO) QCD diagrams as implemented in ARIADNE |]2^ for the QCD 
cascade. As an alternative, samples of events have been generated using the model of lepto 
based on first-order QCD matrix elements plus parton-shower (meps). For the generation of 
the samples with meps, the soft colour interactions option has been switched off. 



In addition, a sample of NC DIS events has been generated using the pythia program |]28 
a lowest-order electroweak calculation including initial- and final-state QCD radiation in the 
leading-logarithm parton-shower approximation. In this case, events have been generated using 
the MRSA pUf set of proton parton densities and the first-order QCD matrix elements have 
not been included. In all cases, the LUND string model ||30| as implemented in jetset |2^ is 
used for modelling the fragmentation into hadrons. 

All MC generated events have been passed through the ZEUS detector and trigger simulation 
programs |T^. They have been reconstructed and analysed by the same program chain as the 
data. 

5 Jet search and energy corrections 

An iterative cone algorithm in the rj — ip plane [jl], is used to reconstruct jets from the 
energy measured in the CAL cells for both data and MC generated events, and also from the 
final-state hadrons for MC generated events. A detailed description of the algorithm can be 



found in |T0|. The jets reconstructed from the CAL cell energies are called cal jets and the 
variables associated with them are denoted by E^l^i, rjl^l and (fill. The axis of the jet is defined 
according to the Snowmass convention 0, where rjl^i (v?ca/) the transverse-energy weighted 
mean pseudorapidity (azimuth) of all the CAL cells belonging to that jet. The energy sharing 
of overlapping jets is dealt with using the following procedure. Two jets are merged if the 
overlapping energy exceeds 75% of the total energy of the jet with the lower energy; otherwise 
two different jets are formed and the common cells are assigned to the nearest jet. The cone 
radius R used in the jet search is set equal to 1. 

For the MC generated events, the same jet algorithm is also applied to the final-state parti- 
cles. The jets found are called hadron jets and the variables associated with them are denoted 
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by E^T%d, VhL and </?^^f^. Hadron jets with E^r^%^ > 14 GeV and -1 < r^fj^ < 2 are selected. 

The comparison of the reconstructed jet variables between the hadron and the cal jets in 
MC generated events |31|] shows no significant systematic shift in the angular variables r^;?^; and 
(flfi with respect to rjl^^^ and ^'had- However, the transverse energy of the cal jet underestimates 
that of the hadron jet by an average amount of ~ 16% with an r.m.s. of 11%. This effect is due 
mainly to energy losses in the inactive material in front of the CAL and is corrected for using the 
following procedure. The transverse energy corrections to cal jets averaged over the azimuthal 
angle are determined using the samples of MC generated events These corrections are 

constructed as multiplicative factors, C*(-ET^ca^ "'^mi); which, when applied to the of the cal 
jets, give the 'corrected' transverse energies of the jets, E^^^ = C{Ej^l^i,T]l^l) x Eij^l^i . 



6 Jet shape 

The differential jet shape is defined as the average fraction of the jet's transverse energy that 
lies inside an annulus in the rj — plane of inner (outer) radius r — Ar/2 (r + Ar/2) concentric 
with the jet defining cone 0: 

(,J: 1 ^ Erjr- Ar/2, r + Ar/2) 

P^'^ N,,,,Arf^^ Et{0,R) ' 

where Exir — Ar /2 , r + Ar /2) is the transverse energy within the given annulus and Njets is the 
total number of jets in the sample. The differential jet shape has been measured for r values 
varying from 0.05 to 0.95 in Ar = 0.1 increments. The integrated jet shape defined by 

is also used. By definition, ip{R) = 1. It has been measured for r values varying from 0.1 to 
1.0 in Ar = 0.1 increments. 

The following procedure is used to reconstruct the differential jet shape from the CAL cells 
in data and MC generated events: for each jet the sum of the transverse energies of the CAL 



cells assigned to the jet, ET,cai{r — Ar/2,r + Ar/2), with a distance r' = ^(Ar/)^ + (Aip)"^ to 
the jet axis between r — Ar/2 and r + Ar/2 is determined and divided by ET,cai{0, 1). The 
differential jet shape as measured with the CAL, Pcaiij), is then defined in analogy with eq. (|lD, 
where the sum now runs over all the cal jets in the selected sample and Njf.ts is the total number 
of cal jets in the sample. Similarly, the integrated jet shape as measured with the CAL, ipcaii^), 
is defined in analogy with eq. (j^). 

The same jet shape definition as used above for the CAL cells is applied to the final-state 
particles in the case of MC generated events and the resulting differential (integrated) jet shape 
is denoted by p^'iir) {ipti^{r))- 



6.1 Jet shape correction 

The differential and integrated jet shapes as measured with the CAL are corrected back to 
the hadron level using the samples of MC generated events. The corrected differential and 
integrated jet shapes, p(r) and ip{r), refer to jets at the hadron level with a cone radius of 
one unit in the rj — (f plane. The measurements are given for jets with Ej^* > 14 GeV and 
— 1 < r]^'^^ < 2 in the kinematic region > 100 GeV^. 
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The corrected jet transverse energy is used only to select the sample of jets {E^r^^ > 14 GeV) 
and to study the dependence of the jet shape as a function of i?jf*. The reconstructed jet shapes 
are then corrected for acceptance and smearing effects using the samples of MC generated 
events. The correction factors also take into account the efficiency of the trigger, the selection 
criteria, the purity and efficiency of the jet reconstruction, and the effects of the energy losses 
due to inactive material in front of the CAL. The corrected differential (integrated) jet shape 
is determined bin-by-bin as p(r) = G^f (r) ■ Pcai{r) and ipi^r) = F^f{r) ■ ipcai{r), where the 
correction factors are defined as G^f (r) = ptdi^)/ pZFir)fnd Fj^f (r) = iptd{r)/^caF{r) 
and are determined separately for each region of r]^'^^ and E'jf . 

For this approach to be valid, the uncorrected jet shapes in the data must be described by 
the MC simulations at the detector level. As shown later, this condition is satisfied by the 
ARIADNE and MEPS simulations in all r]^'^^ and regions studied. The samples of events 
generated with ARIADNE are used to correct the jet shapes. The correction factors G^f{r) do 
not show a strong dependence on t]^^^ or E^?* and vary between 0.7 and 1 for r > 0.15. The 
correction factors for the integrated jet shape F^f{r) differ from unity by less than 25% for 
r > 0.2. Close to the centre of the jet the correction factor G^p(r = 0.05) = F^lf{r = 0.1) is 
large and varies between 1.4 and 1.7 depending on rj^^^ and E^^^. 

The jet shapes have been also reconstructed using tracks instead of CAL cells both in data 
and MC generated events. Since the use of tracks gives an improved spatial resolution for 
the transverse-energy flow of the charged particles within a jet, this study provides a cross- 
check of the resolution in r for the jet shape reconstructed using the CAL. The resulting 
corrected jet shapes are consistent with those using the CAL cells within the uncertainties of 
the measurements (see next section). 



6.2 Systematic uncertainties 

A detailed study of the sources contributing to the systematic uncertainties of the measurements 



has been carried out ||32|. The uncertainties have been classifled into four groups: 



• The energy corrections to the jets and the correction functions to the jet shapes in NC 
and CC DIS have been evaluated using the MEPS generator. The changes induced in p(r) 
are typically below 10%. 

• The absolute energy scale of the cal jets in the MC generated events has been varied by 
±3%. The resulting corrected p(r) changes typically by less than 3%. 

• Variations in the simulation of the CAL response to low-energy particles yielded changes 
in p(r) typically below 3%. 

• Variations in the simulation of the trigger and a variation of the cuts used to select the 
data within the ranges allowed by the comparison between data and MC simulations 
resulted in negligible changes in the corrected jet shapes. 

For the measurements of jet shapes in NC DIS, the statistical errors are negligible compared 
to the systematic uncertainties. Conversely, the statistical errors dominate in the CC DIS 
analysis. The total positive (negative) systematic uncertainty on p(r) at each value of r has 
been determined by adding in quadrature the positive (negative) deviations from the central 
value. The systematic uncertainties have been added in quadrature to the statistical errors and 
are shown as error bars in the figures. 
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7 Results 



7.1 Jet shapes in DIS 

The differential and integrated jet shapes are measured for jets in the reactions 

e+p e+(l7) + jet + X 

with Q'^ > 100 GeV^. Jets are required to have i^^f* > 14 GeV and —1 < t]^^* < 2. There are 
6018, 855 and 53 events in the NC DIS data sample with the range 100 - 1000 GeV^ 

1000 - 5000 GeV^ and 5000 - 25000 GeV^. The corresponding numbers for the CC DIS data 
sample are 84, 123 and 24 events. 

Jet shapes in NC DIS 

The measured differential jet shapes in NC DIS for different regions in r]^^^ and Ej^^ are 
shown in Figures ^ and 0, respectively. The differential jet shape exhibits a prominent peak at 
the centre of the jet. It decreases by a factor ^ 40 from the centre of the jet (r = 0.05) to the 
edge of the jet (r = 0.95). Figure § shows the measured average fraction of the jet's transverse 
energy that lies inside an inner cone of radius r = 0.5 concentric with the jet defining cone, 
ilj{r = 0.5), as functions of t]^'^^ and E^^'^. Note that ilj{r = 0.5) has been measured in ranges of 
Eij^^ and the data points in Figure ^ (lower plot) are located at the weighted mean in each Ej^* 
range. It is observed that the jets become narrower as E^^ increases. The measured 'il>{r = 0.5) 
exhibits no significant dependence on rj^'^^. 

The predictions of ARIADNE, meps and pythia are compared to the measured jet shapes 
in Figures Q to |^. The predicted jet shape of the colour-dipole model (ariadne) describes the 
measured jet shape well in all r/-'^* and E^^ regions considered. The predicted jets of pythia 
tend to be narrower at low E^^^ than those in the data (see Figure |^). In the case of meps, the 
predicted jets show a tendency to be broader at low t]^'^^ than those in the data. 

Jet shapes in CC DIS 

The results for p(r) in CC DIS for different regions of Ej^* are shown in Figure ^. The 
differential jet shape shows similar general features to those of the jets in NC DIS. Figure |] 
shows the measured ^p{r = 0.5) as functions of rj^'^* and Ej^*. The measured ijj{r = 0.5) exhibits 
no significant dependence on rj^^^ and its dependence on E^^'^ is similar to that observed in 
NC DIS. The predictions of ariadne and meps (see Figures ^ and ^) provide a reasonable 
description of the measured jet shape. 

The measured jet shapes in NC DIS are compared with those in CC DIS in Figure |^ and 
found to be very similar in each region of E^^^. Measurements of the ratio of the differential 
jet shapes in CC and NC DIS, p'~^^{r)/p^'^{r), for the same regions of E^* as above are also 
shown in Figure |^ (lower part of each plot) and found to be compatible with unity. In these 
measurements some of the systematic uncertainties common to NC and CC DIS cancel. The 
median of the distribution has been determined for the NC and CC DIS samples of jets 
in each E^'* region: 310 GeV^ (450 GeV^) for 14 < E^^' < 21 GeV, 710 GeV^ (1000 GeV^) 
for 21 < E^"* < 29 GeV, 1260 GeV^ (1600 GeV^) for 29 < E^"* < 37 GeV and 2000 GeV^ 
(2200 GeV^) for 37 < E^^'' < 45 GeV in the NC (CC) DIS samples of jets. Some differences are 
observed in the distributions of the two processes for a given range in ii^jf*. As a cross-check, 
the jet shapes in NC and CC DIS have been measured in a common region of for each range 
in E^* and no significant difference has been found. Therefore, the observation that the jet 
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shapes in NC and CC DIS are very similar, for the same range of -E"^^*, is independent of the 
different distributions in these processes. 

7.2 Comparison to jet shapes in photoproduction 

In photoproduction, two types of QCD processes contribute to jet production at LO |3^ : 
either the photon interacts directly with a parton in the proton (direct process) or the photon 
acts as a source of partons which interact with those in the proton (resolved process). It has been 
noted that resolved processes dominate jet photoproduction in the entire E^r^^ region studied 
3T| . In the case of dijet photoproduction the contributions of resolved and direct processes can 



be separated |^ by using the variable x'^^^ = (E.eis ^^'^e"^''*) 7(2^7), where the sum runs 



J jets ■ 

over the two jets of highest E^^^ and E^ is the initial photon energy. This variable represents 
the fraction of the photon's momentum participating in the production of the two jets with 
highest Ej^^. The LO direct and resolved processes largely populate different regions of x^^^ , 
with the direct processes being concentrated at high values. 

In Figure |^ the measured integrated jet shape in NC DIS are compared to those in dijet 
photoproduction ||T^ for two different regions: x^^^ > 0.75 and x^^^ < 0.75. The comparison 
between the jet shapes in NC DIS and dijet photoproduction is made for the same ranges of 
77-'^* and the Ej^* spectrum is similar in these two processes. The jets produced in NC DIS are 
narrower than those in dijet photoproduction but closer to those dominated by direct processes 
^^oss ^ 0.75). This comparison can be understood in terms of the large fraction of final- 
state quark jets expected in NC DIS {e~^q e~^q) and direct processes in photoproduction 
(dominated by the subprocess 'jg — >■ qq). The remaining differences may be attributed to the 
contribution from the direct subprocess jq — > qg and that of resolved processes, in which the 
jets are broader as shown by the measurements in dijet photoproduction with x^^^ < 0.75. 

7.3 Comparison to measurements in e^e~ and pp collisions 

The measured jet shape in NC (CC) DIS with > 100 GeV^ for jets with transverse energy 
between 37 and 45 GeV, with a mean of 40 GeV (41 GeV), is compared to the measurements 
of the jet shape corrected to the hadron level in pp collisions by CDF ^ and D0 [0 and in 
e"'"e~ interactions by OPAL 0: 

• The CDF data have been obtained using an iterative cone algorithm with R = 1 
similar to that used here. The measurements shown are for jets with transverse energy 
between 40 and 60 GeV, with a mean of 45 GeV, and pseudorapidity 0.1 < {rj^'^^l < 0.7. 
The contribution to the jet shape due to the underlying event was found to be small^. If 
a jet shares more than 75% of its energy with a jet of higher energy, the two are merged 
together; otherwise, they are defined as distinct and the particles common to both jets 
are assigned to the nearest jet. 

• The D0 data |^ have been obtained also using an iterative cone algorithm with R = 1 
similar to that used here. The jet direction was defined according to a convention different 
from that of Snowmass; however, this difference is not expected to have a significant 



effect on the results ||36|. The jet shape has been measured for jets with transverse energy 
between 45 and 70 GeV, with a mean of 53 GeV, and pseudorapidity \7]^^^\ < 0.2. The jet 
shape has been corrected to remove the small contribution due to the underlying event. 



^In NC and CC DIS the underlying event is not expected to contribute in the kinematic region studied here. 
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Two jets were merged if more than 50% of the Et of the jet with smaller was contained 
in the overlap region; otherwise, the two jets were not merged and each particle in the 
overlap region was assigned to the nearest jet. 

• The OPAL data have been obtained using a cone algorithm especially designed to 
emulate that of the CDF measurements, i.e. defining the cone in the f] — <f plane, using 
R = 1, demanding \ri^'^^\ < 0.7 and measuring the transverse energy flow. The jet shape 
has been measured for jets with energy greater than 35 GeV, with a mean of 40.4 GeV. 
The e+e~ data have no underlying event. Overlapping jets are treated using the same 
procedure as CDF. 

The measured differential jet shapes in NC and CC DIS are compared to that measured 
in e^e~ interactions in Figure ^ and are found to be similar. The ratio of the differential jet 
shapes in NC DIS and e~^e~ interactions, p^'^(r)/p^^^ (r), is also shown in Figure § (lower 
part of the figure) and is found to be compatible with unity within the uncertainties of the DIS 
measurements, which are dominant. For the selected samples of jets, the jet shapes in e^e^ 
interactions and DIS are expected to be similar due to the large fraction of final-state quark jets 
in these two processes. However, some differences may appear since there are configurations of 
colour flow (for example, that of initial-state QCD radiation) in DIS which are not present in 
e~^e~ . The striking similarity in the jet shapes indicates the large extent to which the pattern 
of QCD radiation within a quark jet is independent of the hard scattering process in these 
reactions. 

The measured integrated jet shapes in DIS are compared to those in e+e~ interactions and pp 
collisions in Figure ^ The measured jets in DIS at HERA are found to be narrower than those 
in pp collisions. The measurements in pp collisions have been performed for jets with slightly 
higher energy than those in NC and CC DIS. This difference cannot explain the discrepancy 
in the jet shapes since the jets become narrower as the jet energy increases. As stated in 
most of the difference between the jet shapes in e~^e~ interactions and pp collisions can be 
ascribed to the larger fraction of gluon jets in the latter reaction. The comparison between the 
measured jet shapes in DIS and pp collisions suggests that, also in this case, the difference can 
be attributed to differences between quark and gluon jet properties. 

8 Summary and conclusions 

Measurements have been presented of the differential and integrated jet shapes in neutral- and 
charged- current deep inelastic e^p scattering at -^i = 300 GeV using data collected by ZEUS 
in 1995 and 1996. The jet shapes refer to jets at the hadron level with a cone radius of one 
unit in the Tj — if plane and are given for the kinematic region > 100 GeV^. Jets with 
Ej^^ > 14 GeV and —1 < rj^^* < 2 have been considered. The jets become narrower as 
increases. No significant rj^'^^ dependence of the jet shape has been observed. The measured jet 
shapes in neutral- and charged-current DIS are found to be very similar. 

The measurements of jet shapes have been compared to the predictions of Monte Carlo gen- 
erators using different models for the QCD radiation. The colour-dipole model as implemented 
in ARIADNE provides a reasonable description of the measured jet shapes in all t]^^^ and E^* 
regions studied. The parton-shower approach without first-order QCD matrix-elements predict 
jets which are slightly narrower at low -Ejf* than those in the data for all the rj^^* regions stud- 
ied. The inclusion of first-order QCD matrix-elements improves the description of the data for 
-qjet ^ leads to jets which are slightly broader for rj^'^'^ < 1. 
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The jets in neutral-current DIS are narrower than those in dijet photoproduction but closer to 

those in direct-photon processes for the same ranges in jet transverse energy and pseudorapidity. 
The jets in DIS are found to be narrower than those in pp collisions. This difference can be 
attributed to a larger contribution of gluon jets in pp collisions. The measured jet shapes in 
neutral- and charged-current DIS are similar to those in e'^e" interactions for comparable ranges 
of jet transverse energy. Since the jets in e+e" interactions and deep inelastic e+p scattering 
are predominantly quark initiated, the similarity in the jet shapes indicates that the pattern 
of QCD radiation within a quark jet is to a large extent independent of the hard scattering 
process in these reactions. 
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Figure 1: Measured differential jet shapes corrected to the hadron level, p(r), in neutral-current 
DIS with > 100 GeV^ for jets with E^^^ above 14 GeV in different r/^"* regions (black 
dots). The error bars include the statistical and systematic uncertainties added in quadrature 
(typically smaller than the dots). The predictions of PYTHIA (dotted lines), ARIADNE (solid 
lines), and MEPS (dashed lines) are shown for comparison. The predictions have been obtained 
by an integration over the same bins as for the data and are presented as smooth curves joining 
the calculated points. 
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Figure 2: Measured differential jet shapes corrected to the hadron level, p(r), in neutral-current 
DIS with > 100 GeV^ for jets in the rj^^'^ range between —1 and 2 in different Ei^'^ regions 
(black dots). The error bars include the statistical and systematic errors added in quadrature. 
The predictions of pythia (dotted hnes), ARIADNE (sohd lines), and MEPS (dashed hnes) are 
shown for comparison. The predictions have been obtained by an integration over the same 
bins as for the data and are presented as smooth curves joining the calculated points. 
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Figure 3: The measured integrated jet shape corrected to the hadron level at a fixed value of 
r = 0.5, ijj{r = 0.5), as a function of 77^^* (upper plot) and E^* (lower plot), in neutral-current 
DIS with g2 > 100 GeV^ for jets with S^^* > 14 GeV in the r^^'^* range between -1 and 2 
(black dots). The error bars include the statistical and systematic errors added in quadrature. 
The predictions of PYTHIA (dotted lines), ARIADNE (solid lines), and MEPS (dashed lines) are 
shown for comparison. The predictions have been obtained by an integration over the same 
bins as for the data and are presented as smooth curves joining the calculated points. 
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Figure 4: Measured differential jet sliapes corrected to tfie fiadron level, p(r), in charged-current 
DIS with > 100 GeV^ for jets with -1 < ry^'^* < 2 in different * regions (black dots). The 

predictions of ARIADNE (solid lines) and MEPS (dashed lines) arc shown for comparison. The 
predictions have been obtained by an integration over the same bins as for the data and are 
presented as smooth curves joining the calculated points. 
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Figure 5: The measured integrated jet shape corrected to the hadron level at a fixed value of 
r = 0.5, 'ilj{r = 0.5), as a function of 77-^*^* (upper plot) and -E^^* (lower plot), in charged- current 
DIS with Q2 > 100 GeV^ for jets with E^^^^ > 14 GeV in the r]^^* range between -1 and 2 (black 
dots). The inner error bars represent the statistical errors of the data, and the outer errors 
bars show the statistical and systematic uncertainties added in quadrature. The predictions 
of ARIADNE (solid lines) and meps (dashed fines) are shown for comparison. The predictions 
have been obtained by an integration over the same bins as for the data and are presented as 
smooth curves joining the calculated points. 
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Figure 6: Measured differential jet shapes corrected to the hadron level, p(r), in charged- current 
DIS with g2 > 100 GeV^ for jets with -1 < 77^'^* < 2 in different E^j.^^ regions (open circles). 
The measured jet shapes corrected to the hadron level for jets in neutral-current DIS with 

> 100 GeV^ with — 1 < t]-'^* < 2 are shown for comparison (open squares). Measurements 
of the ratio p'^'^ (r) / p'^^ (r) arc shown underneath each plot. The inner error bars represent 
the statistical errors of the data, and the outer errors bars show the statistical and systematic 
uncertainties added in quadrature. 
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Figure 7: Measured integrated jet shape corrected to the hadron level, in neutral-current 
DIS with Q2 > 100 GeV^ for jets with above 14 GeV and -1 < r]^^^ < 2 (black dots). 
The measured jet shape corrected to the hadron level for jets in dijet photoproduction with 



Ej?'^ above 



14 GeV and — 1 < rj^ < 2 is shown for comparison: for dijet production with 



X, 



OBS 



< 0.75 (stars) and for dijet production with > 0.75 (open circles). 
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Figure 8: Measured differential jet shapes corrected to the hadron level, p{r), in neutral- 
(charged-) current DIS with > iqo GcV^ and a median of 2000 GeV^ (2200 GeV^) for jets 
with T]^'^* in the range between —1 and 2 and 37 < E^^ < 45 GeV are shown as squares (open 
circles). The measurements in CC DIS have been obtained for the same values of r as those in 
NC DIS, and for an easier comparison the measurements are plotted at r + 0.025. The mea- 
surements of the jet shape in e'^e~ interactions by OPAL (black dots) is shown for comparison. 
The ratio of differential jet shapes in NC DIS and e'^e~ interactions, p^^ {r)/ p^^^ (r), is shown 
in the lower part of the figure. 
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Figure 9: Measured integrated jet shapes corrected to the hadron level, iplr), in neutral- 
(charged-) current DIS with Q'^ > 100 GeV^ for jets with rj^^* in the range between —1 and 2 

and 37 < E^^^^ < 45 GeV are shown as stars (diamonds). The measurements of jet shapes in 
pp collisions by CDF (triangles) and D0 (circles) and in e~^e~ interactions by OPAL (squares) 
are shown for comparison. 
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